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Abstract 

We investigate how the CP quantum numbers of a neutral Higgs boson or spin-zero resonance 
<I>, produced at the CERN Large Hadron Collider, can be determined in its T-pair decay mode 
4> — > t~t + . We use a method [1] based on the distributions of two angles and apply it to the 
major 1 -prong x decays. We show for the resulting dilepton, lepton-pion, and two-pion final 
states that appropriate selection cuts significantly enhance the discriminating power of these 
observables. From our analysis we conclude that, provided a Higgs boson will be found at 
the LHC, it appears feasible to collect the event numbers needed to discriminate between a 
CP-even and CP-odd Higgs boson and/or between Higgs boson(s) with CP-conserving and 
CP-violating couplings after several years of high-luminosity runs. 
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I. INTRODUCTION 



If the Large Hadron Collider (LHC) at CERN will reach its major physics goal of discov- 
ering a spin-zero resonance, the next step will be to clarify the question whether this is the 
standard model (SM) Higgs-boson or some nonstandard resonance, as predicted by many of 
the presently discussed new physics scenarios. (For reviews, see 10-0] •) Part of the answer to 
this question will be given by measuring the CP quantum numbers of such a particle. There 
have been a number of proposals and investigations on how to determine these quantum 
numbers for Higgs-like resonances 4>, for several production and decay processes at hadron 



colliders, including Refs. IUU7H19H. (For an overview, see |20D.) It is the purpose of this ar- 
ticle to study a method |l]] with which one can pin down whether such a state is CP-even, 
CP-odd, or a CP- mixture, namely in its decays into T-lepton pairs with subsequent 1 -prong 
decays. 

Our investigations below are applicable to neutral spin-zero resonances hj, for instance to the 
Higgs-boson(s) of the standard model and extensions thereof, with flavor-diagonal couplings 
to quarks and leptons as described by the Yukawa Lagrangian 

J? Y = -(V2G F ) 1 / 2 Y,mf(a,fff + b jf fiy 5 f)h j . (1) 

j,f 

Here is the mass of the fermion / and we normalize the coupling constants to the Fermi 

constant Gp. A specific model is selected by prescribing the reduced scalar and pseudoscalar 
Yukawa coupling constants ajf and bjf. In the standard model (SM) with its sole Higgs- 
boson, 7=1 and ajf = 1, bjf = 0. Many SM extensions predict more than one neutral spin- 
zero state and the couplings (OQ) can have more general values. Two-Higgs doublet models, 
for instance, the nonsupersymmetric type-II models and the minimal supersymmetric SM 
extension (MSSM, see, e.g., [Q-Q, 01) contain three physical neutral Higgs fields hj. If 
the Higgs sector of these models is CP-conserving, or if Higgs-sector CP violation (CPV) 
is negligibly small, then the fields hj describe two scalar states, usually denoted by h and 
H, with bjf = 0, ajf ^ 0, and a pseudoscalar, denoted by A, with ajf = 0, bjf ^ 0. In 
the case of Higgs-sector CPV, the mass eigenstates hj are CP mixtures and have nonzero 
couplings ajf ^ and bjf ^ to scalar and pseudoscalar fermion currents. This would lead 
to CP- violating effects in the decays hj — >■ // already at Born level [0]. 
In the following, we use the generic symbol $ for any of the neutral Higgs-bosons hj of the 
models mentioned above or, in more general terms, for a neutral spin-zero resonance. At 
the LHC, a resonance can be produced, for instance, in the gluon and gauge boson fusion 
processes gg — > <& and qiqj — > Qq\q ! j, as well as in association with a heavy quark pair, tt<3> 
or bb$>. Recent studies on Higgs-boson production and decay into x leptons within the SM 



and the MSSM include 112 11 12211 . Our method for determining the CP parity of 4> can be 
applied to these and to any other LHC ^-production processes. 

The spin of a resonance 4> can be inferred from the polar angle distribution of the 4>-decay 
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products. In its decays to T leptons, <E> — > t~t + , which is apromising LHC search channel 
for a number of nonstandard Higgs scenarios (see, e.g. [Q, @] and the recent LHC searches 
I23 , 24]) T-spin correlations induce specific angular distributions and correlations between 
the directions of flight of the charged T-decay products, in particular an opening angle distri- 
bution and a CP-odd triple correlation and associated asymmetries Once a resonance 
is discovered, these observables can be used to determine whether it is a scalar, a pseu- 
doscalar, or a CP mixture. 

The discriminating power of these observables can be exploited fully if the x ± rest frames, 
i.e., the x energies and three-momenta can be reconstructed. At the LHC this is possible for x 
decays into three charged-pions. With these decay modes the CP properties of a Higgs-boson 
resonance can be pinned down efficiently II 1811 . For x decays into one charged particle the 
determination of the x ± rest frames is, in general, not possible at the LHC. For the 1-pron 
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decays x ± —t a ± the a+aT zero-momentum frame can, however, be reconstructed. In 
two observables, to be determined in this frame, were proposed and it was shown, for the 
direct decays T + T~ — > n + n~v T v T , that the joint measurement of these two observables de- 
termines the CP nature of It will even be possible to distinguish (nearly) mass -degenerate 
scalar and pseudoscalar Higgs-bosons with CP-invariant couplings from one or several CP 
mixtures. 

In this paper we analyze the T-pair decay mode of <I> for all major 1 -prong T-decay s x ± — > a ± 
and investigate how this significantly larger sample can be used for pinning down the CP 
quantum numbers of <J> in an efficient way. As the respective observables originate from 
T-spin correlations, the T-spin analyzing power of the charged particle a is crucial for this 
determination. The T-spin analyzing power of the charged lepton in the leptonic T decays 
and of the charged-pion in the 1 -prong hadronic decays x ± — > p ± ,af — > is rather poor 
when integrated over the energy spectrum of the respective charged prong. Thus, the crucial 
question in this context is whether experimentally realizable cuts can be found such that, on 
the one hand, the T-spin analyzing power of the charged prongs is significantly enhanced 
and, on the other hand, the data sample is not severely reduced by these cuts. We have 
studied this question in detail and found a positive answer. 

The paper is organized as follows. In the next section we briefly describe the matrix elements 
on which our Monte Carlo event simulation is based, and we recapitulate the two observables 
with which the CP nature of a Higgs-boson can be unraveled. In Sec. [Ill] we analyze in detail 
the distribution that discriminates between a scalar and pseudoscalar Higgs-boson, both for 
dilepton, lepton-pion, and two-pion final states, for several cuts. We demonstrate for a set 
of "realistic", i.e., experimentally realizable cuts that the objectives formulated above can 
actually be met. This is then also shown for the distribution that discriminates between 4> 
bosons with CP- violating and -conserving couplings. We conclude in Sec.lTVl 
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II. DIFFERENTIAL CROSS SECTION AND OBSERVABLES 



We consider the production of a spin-zero resonance $ - in the following collectively called 
a Higgs-boson - at the LHC, and its decay to a pair of T ± leptons: 

pp^<f> + X^T~T + +X. (2) 

The decays of x ± are dominated by 1 -prong modes with an electron, muon, or charged-pion 
in the final state. We take into account the following modes, which comprise the majority of 
the 1 -prong x decays: 



T -> Z + V/ + V 

t a x + v T n + 2n° + v T , 

z n + v T . (3) 

In the following = / T , 7T T refer to the charged prongs in the decays ©. 
The hadronic differential cross section do for the combined production and decay processes 
©, © can be written as a convolution of parton distribution functions and the partonic 
differential cross section do for P1P2 — > 4> — )• t + t~ — > a + a'~ + X (where p\ and p2 are 
gluons or (anti)quarks): 

x I A + b(E a ,-)B + -q--b(E a+ )B--q+-b(E a ,-)b(E a+ ) £ C (7 ^+ 



Here, y/s is the partonic center-of-mass energy, f5 r = \J I — \m\jp\, 

Y,\M(piP2^<Z> + X)\ 2 and D~ l (<*>) = (p| -m| + im$r^) _1 (5) 

is the spin and color averaged squared production matrix element and the Higgs-boson prop- 
agator, respectively, with m<j>, p^ and rj£ £ denoting the Higgs-boson mass, its 4-momentum 
and its total width 4 . The squared matrix element \T\ 2 of the decay — > t + t~X, integrated 
overX, is of the form 

\T\ 2 = V2G F m 2 T (A+B+sf + BTsr+QjS+sr) , (6) 



4 Equation © holds as long as nonfactorizable radiative corrections that connect the production and decay 
stage of <t> are neglected. 
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Table I: Tree-level coefficients of the squared decay matrix element ©, (Q for <J> = H, A (scalar, 
pseudoscalar) and for a CP mixture. 



where are the normalized x ± spin vectors in the respective x ± rest frames. The dynamics 



± 



of the decay is encoded in the coefficients A, Bf and Qj. Rotational invariance implies that 



B 



± 



Qj = C\8ij + C2&,- + C3 



/ ' 



(7) 



where k (k ) is the (normalized) x momentum in the t + t zero-momentum frame 
(ZMF). At tree level, B ± = 0. (A nonzero absorptive part of the amplitude, induced for 
instance by the photonic corrections to $ — > xx renders these coefficients nonzero, but the 
effect is very small lll4ll .) The tree-level coefficients A and ci,2,3 induced by the general 
Yukawa couplings © are given in Table [T] (cf . also lll4ll ) for a scalar (b t = 0) and a pseu- 
doscalar (a T = 0) Higgs-boson, 4> = H, A, and a CP mixture. 

We use the narrow-width approximation for T*. The branching ratios of the 1 -prong x decays 
© are denoted by Br ± ± = T x ±^ a ± /T 1 " 1 . Moreover, the measure dQ. x = d cos Q x d<px in © 
is the differential solid angle of the T~ in the Higgs rest frame, and d£l a ± = dcos6 a ±d(p a ± 
and E a ± is the differential solid angle and the energy of the charged prong a ± in the x ± rest 
frame. Furthermore, the functions n(E a ^) and b(E a ^) encode the decay spectrum of the 
respective 1 -prong polarized T T decay and are defined in the T T rest frames by 



dT(T T (s T ) ->a T (^)+X) 
T(tT ->■ +X) dE a ^dQ. a ^ / '(Alt) 



n{E a ^(\±b(E^)^-iF) 



(8) 



where q T is the normalized momentum vector of the charged prong in the respective 
frame. The function n{E a ) determines the decay rate of x — > a while b(E a ) encodes the 
T-spin analyzing power of the charged prong a = /, 7T. We call them spectral functions for 
short. They are given for the T-decay modes © in Appendix A. Using the spin-density 
matrix formalism, the combination of © and © yields, with ©, the formula ©. 
The decay distribution © and the coefficients c\, C3 of TableUimply that, at the level of the 
T + T~ intermediates states, the spin observables s + ■ s~ and k ■ (s + x s~ ) discriminate between 
a CP-even and CP-odd Higgs-boson, and between a Higgs-boson with CP-conserving and 
CP- violating couplings, respectively [@,[l4]. A the level of the charged prongs a + a'~, these 
correlations induce a nontrivial distribution of the opening angle Z(q + , q~) and the CP-odd 
triple correlation k ■ (q + x q~), as can be read off from ©. The strength of these correlations 
depends on the product b(E a i-)b(E a +), while n(E a r-)n(E a +) is jointly responsible for the 
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number of a + a'~ events 5 . 

A direct analysis of experimental data in terms of the kinematic variables used in the dif- 
ferential cross section Eq. © is not possible since the momenta of the x decay products are 
measured in the laboratory frame and the reconstruction of the x ± and <I> rest frames is, in 
general, not possible. In Ref. yj] it was shown that one can, nevertheless, construct exper- 
imentally accessible observables with a high sensitivity to the CP quantum numbers of 4>. 
The crucial point is to employ the zero-momentum frame of the a + a'~ pair. 
The distribution of the angle 

<p* = arccos(n^_ + -n^) (9) 

discriminates between a J PC = ++ and h state. Here h*^ are normalized impact param- 
eter vectors defined in the zero-momentum frame of the a + a'~ pair. These vectors can be 
reconstructed 111] from the impact parameter vectors n T measured in the laboratory frame by 
boosting the 4- vectors = (0,n=p) into the a'~a + ZMF and decomposing the spatial part 
of the resulting 4-vectors into their components parallel and perpendicular to the respective 
7T T or / T momentum. We emphasize that (p* defined in Eq. © is not the true angle between 
the x decay planes, but nevertheless, it carries enough information to discriminate between a 
CP-even and CP-odd Higgs-boson. 

The role of the CP-odd and P-odd triple correlation mentioned above is taken over by the 
triple correlation Gq P = p*l ■ (n^ + x n^~) between the impact parameter vectors just defined 
and the normalized a'~ momentum in the a'~a + ZMF, which is denoted by p!_. Equivalently, 
one can determine the distribution of the angle |l]] 

y/£ P = arccos(p* • (n^+ x A*-)) . (10) 

In an ideal experiment, where the energies of the x decay products a in the x ± rest frames 
would be known, one could determine the coefficients A, B ± , and ci 2,3 by fitting the dif- 
ferential distribution © (using the SM input of the Appendix) to the data. However, due to 
missing energy in the final state, detector resolution effects and limited statistics, one has to 
average over energy bins. Moreover, for ^ 7T T the function b(E) is not positive (nega- 
tive) definite, see below. Therefore, energy averaging can lead to a strong reduction of the 
sensitivity to the coefficients of b(E) in the differential cross section. A judicious choice of 
bins or cuts is therefore crucial to obtain maximal information on the CP properties of 4>. 
We will discuss this in detail in the next section. 



5 The integral J n(E a )b{E a )dE a determines the overall T-spin analyzing-power of the particle a. It seems 
worth recalling that the physics of T decays, i.e., the V — A law, has been tested to a level of precision which 
is much higher than what is needed for our purposes. Therefore, when comparing predictions with future 
data, one can use the functions n(E a ) and b(E a ) of the Appendix as determined within the standard model. 
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III. RESULTS 



The observables © and (flOl) can be used for the 1 -prong T-pair decay channels of any Higgs- 
boson production process at the LHC. We are interested here in the normalized distributions 
of these variables. If no detector cuts are applied, these distributions do not depend on the 
momentum of the Higgs-boson in the laboratory frame; i.e., these distributions are indepen- 
dent of the specific Higgs-boson production mode. Applying selection cuts, we have checked 
for some production modes (see below) that, for a given Higgs-boson mass m<6 > 120 GeV, 
the normalized distributions remain essentially process-independent (see also [ lj]). 
For definiteness, we consider in the following the production of one spin-zero resonance $ at 
the LHC (VS= 14TeV) in a range of masses m<£ between 120 and 400 GeV. As we employ 
the general Yukawa couplings (OQ), our analysis below can be applied to a large class of 
models, including the standard model, type-II 2-Higgs doublet models, and the Higgs sector 
of the MSSM. Within a wide parameter range of type-II models, 4> production is dominated 
by gluon-gluon fusion; for large values of the parameter tan/3 = V2/V1 (where vi 5 2 are the 
vacuum expectation values of the two Higgs doublet fields) the reaction bb — > $ takes over. 
(For a recent overview of various Higgs-boson production processes and the state-of-the-art 
of the theoretical predictions, see, e.g., 11251 26]. Higgs-boson production and decay into 
T~T + was analyzed in the SM and MSSM in [21,[22|], taking recent experimental constraints 
into account.) 

For obtaining the results given below we have used the production processes bb — >■ 4> and 
gg — > <I>. The reaction chains ©, © were computed using leading-order matrix elements 
only, but our conclusions will not change when higher-order QCD corrections are taken into 
account or other Higgs-boson production channels with large transverse momentum pf are 
considered. Our method can be applied to all production channels, because no reconstruction 
of the Higgs-boson momentum or the x momenta is needed for the determination of the 
distributions © and (PTOl) . Therefore, the method is applicable to Higgs-boson production 
with small or large transverse momentum, as long as the Higgs resonance can be identified 
in the %% events. (For a discussion of the background see the end of Sec. IIII DP 
If pj is small, the distributions can be measured as described below. If the v~, T + decay into 
leptons or via a p or a\ meson, an approximate reconstruction of the Higgs-boson rest frame, 
as outlined in Sec. IIII C[ will increase the discriminating power of the distributions, because 
appropriate cuts in this frame separate T-decay particles with large and small energies. 
If pj is large, the reconstruction of the Higgs rest frame can be performed, see |27j]. With 
similar cuts as used below, this leads to an even better discriminating power of the (p* and 
distributions. 

We have implemented Eq. © into a Monte Carlo simulation program which allows us to 
study the reconstruction of observables in a variety of reference frames and to impose selec- 
tion cuts on momenta and energies. 

In Sec. IIII Al - IIII Dl we analyze, for the various 1-prong final states, the (p* distributions for 
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(a) T -> 1 + V, + v T (b) $ -> T ~T + -» + 3v 




Ei 



Figure 1: (a) The spectral functions n(E/) and b{E{), Eq. (I22l ). for the leptonic x decay. The function 
n(Ei) is given in units of GeV -1 . (b) The normalized <p* distribution for In final states without 
selection cuts in the laboratory frame for a Higgs mass of m<& = 400 GeV. A cut on the lepton energy 
in the X rest frame at m T /4 ~ 0.44 GeV serves to show the effect of rejecting events where b{E{) is 
positive and negative, respectively. 

a scalar and pseudoscalar Higgs-boson, i.e., a spin-zero resonance $ with reduced Yukawa 
couplings a T ^ 0,Z? T = and a x = 0,b T ^ 0, respectively, to X leptons. For definiteness we 
choose a T = 1 and b r = 1, respectively. The distribution of the CP angle ty/£p is computed in 
Sec. IIII El for Higgs-bosons with CP-violating and CP-conserving couplings. 

A. Lepton-pion final state: xx — > 1% + 3v 

We start by discussing the case where the x~ from — > x~x + decays leptonically, x~ — > 
V + V/ + v T , and the x + undergoes a direct decay into a pion, x + — > K + + v T . The purpose of 
this section is to study the shapes of the (p* distributions for scalar and pseudoscalar Higgs- 
bosons when cuts are applied to the charged lepton; therefore, no cuts are applied at this 
point to the pion energy and momentum. 

The charged lepton energy spectrum in the x — >■ / decay is determined by the functions n(E[) 
and b(Ej) given in the Appendix, Eq. (1221) . These functions are shown in Fig. a). One 
sees that the function b(Ei), which determines the T-spin analyzing power of /, changes sign 
at E[ = m T /4. Therefore also the slope of the resulting (p* distribution for nl final states 
changes sign at this energy. The optimal way to separate a CP-even and CP-odd Higgs- 
boson would be to separately integrate over the energy ranges E\ > m T /4 and E\ < m T /4. 
The resulting (p* distributions for a scalar (H, red lines 6 ) and pseudoscalar (A, black lines) 

6 Color in the electronic version. 
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Figure 2: Normalized lepton energy distribution (in the T rest frame) for different Higgs-boson 
masses, with and without selection cuts fTTV 

Higgs-boson are shown in Fig.[TJb). For the energy range < E\ < m T /4, the <p* distribution 
has a positive (negative) slope for a scalar (pseudoscalar) Higgs-boson (dashed curves). For 
Ei > m T /4 the slopes change sign and the difference between a CP-even and a CP-odd boson 
becomes more pronounced (solid curves). However, at a LHC experiment, the separation of 
these two energy ranges is not possible because the x momenta can not be reconstructed and, 
therefore, the lepton energy E\ in the x rest frame can not be determined. 
The difference between the (p* distributions for a scalar and pseudoscalar Higgs-boson is, 
however, not completely washed out by integrating over the full lepton energy range, be- 
cause both n(E[) and b{E{) have a significant energy dependence, see Fig. [Da). The region 
< Ei < m T /4 contributes only about 19% to the decay rate r T ^/; in addition, b{E{) is 
small in this energy range. Therefore, after integration over the full £/ range, the <p* distri- 
bution is already quite close to the solid lines of Fig. \V[b). Moreover, one can suppress the 
contribution from the low-energy part of the spectrum by imposing a cut on the transverse 
momentum of the charged lepton in the laboratory frame. For the LHC experiments, suitable 
selection cuts on the transverse momentum and the pseudorapidity of the lepton in the pp 



The effect of these cuts on the normalized lepton energy distribution in the x~ rest frame is 
shown in Fig. [21 Rejecting events with small p l T preferentially removes events with small 
lepton energy in the X rest frame. The effect is more pronounced for light Higgs-boson 
masses where the x energy is smaller on average. For m$ = 120 GeV, only a small fraction 
of ill events with E { < m T /4, about 3.6%, survives the cuts (ITTb . For m$ = 200 and 400 GeV 
the corresponding fractions are 9.4% and 14%, respectively. Events with Ei < m t /4 that 



frame are |23j,|24j]: 




(11) 
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Figure 3: The normalized <p* distributions for Ik3v final states. The solid (dashed) curves show the 
distribution with the cuts p l T > 20 GeV and \r\i\ < 2.5 (without these cuts), (a) m$ = 120 GeV; (b) 
mo = 400 GeV. 

pass the above cuts have energies close to m z /4, where the function b(E[) is very small. The 
resulting <p* distribution is almost unaffected by contributions with E\ < m % /4, for Higgs 
masses up to 200 GeV. As an example, the <p* distributions are displayed for m$ = 120 GeV 
in Fig.Oa) and for = 400 GeV in Fig.[3£b). From these results we conclude that only for 
very large Higgs-boson masses one can expect to improve the discrimination of scalar and 
pseudoscalar bosons by such a detector cut. 

The experimentally relevant case, where in addition also selection cuts on the charged-pion 
are applied, will be discussed in Section HITDl 



B. Hadronic final states: t r + ^-{a l ,p ,71 }n 



Next we analyze the case where the x~ decays to %~ either via a p meson, T~ — > p + v T — > 
7l~ + 7T° + V T , an a\ meson, T~ — y ci^ -\- v-r — y K~ -\- 27t° + v T , or directly, -> K~ + v T , 
while T + undergoes a direct 2-body decay, T + — > TZ + + v T . (The respective branching ratios 
are collected in Table PVT). 

The spectral functions n{E n ) and b(E n ) for the p and a\ modes are given in the Appendix 
and shown in Fig. @J The direct decay mode T~ — > %~ + v T is characterized by a constant 
pion energy in the x rest frame and has maximal T-spin analyzing power b = 1 . 
As in the case of leptonic x decay the functions b p (E K ) and b ai (E n ) change sign, at approx- 
imately 0.55 GeV. In contrast to the leptonic case, however, contributions from small pion 
energies are not suppressed by small differential rates, as evidenced by the functions n p (E 7t ) 
and n ai (E K ). At the LHC one will probably not be able to distinguish between the different 
1 -prong decay modes into a pion, at least not in an efficient way. Thus, one has to combine in 
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Figure 4: Charged-pion spectral functions n{E n ) and b(E % ) for hadronic X decays: (a) X — >• p V T — >• 
7T~7T v T ; (b) T~ — >• a7v T — > 7r _ 27r°v T . The functions and n{E n )b{E n ) are given in units of 

GeV" 1 . 

(a) t" -> (ar.pl + v T -> tT + X (b) * -> t"t + -» (af.p - )^ + 2v -> ;T7r + + X 




Figure 5: (a) Spectral function n(E J[ ) for the combined hadronic decay modes X — > a\ ,p — > 7t (solid), 
(b) The distribution of the charged pion energy E n for the combined a[ and p ~ decays, for different 
Higgs-boson masses, with and without selection cuts. 



the Monte Carlo modeling the different decay modes, weighted with their branching ratios. 
The combined functions n(E K ) for the decays T — > a\ and x — > p are shown in Fig. Oa). 
The combined distribution is dominated by the p decay mode because of its larger branching 
ratio. The result (solid red curve) in Fig.[5|a) shows that contributions where E K > 0.55 GeV 
and E n < 0.55 GeV, i.e., where b{E n ) is positive and negative, respectively, are equally 
important. Without cuts, one would, as a consequence, not be able to distinguish between 
scalar and pseudoscalar Higgs-bosons by means of the (p* distribution. Therefore we im- 
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Figure 6: The normalized <p* distributions for the combined decays <I> — > T~T + — > {a^ ,p~}7i + — > 
7t~7t + . The px and T] cuts (IT2T > are imposed on 7T~ only (dashed curves). The solid curves show the 
distributions which result if instead of (PT2l an experimentally not feasible cut on the pion energy E n - 
in the T rest frame is applied, (a) m$ = 200 GeV; (b) m$ = 400 GeV. 

pose the following selection cuts on the charged-pion in the pp laboratory frame, which are 



compatible 7 with the cuts used by the LHC experiments II23L 12411 



/>£>40GeV, |%|<2.5. (12) 

At this point, the constraints (fl~2l) are imposed - for the purpose of analyzing the effect of 
these cuts - only on the %~ from T~ decays. These cuts will be applied to both %~ and n + 
from T~ and T + decays, respectively, in Section ITODl 

The impact of these cuts on the distribution o [ da /dE n (where E n is the energy of the 7T~ 
in the T~ rest frame) for the combined T~ — > a^,p^ — > n~ decay modes is displayed in 
Fig. Ob) for several Higgs-boson masses between 120 and 400 GeV. The solid curve shows 
the distribution without cuts; it does not depend on m<$. The cuts (fT2l) preferably reject 
events with small E K . The effect of these cuts is strong for light Higgs-boson masses and 
still pronounced for m<j ~ 200 GeV. 

As a consequence, the (p* distributions are dominated by contributions with positive values of 
the T-spin analyzer functions b p (E n ), b ai (E n ), and the distributions clearly differ for ^> = H 
and = A, see Fig.[6ta). For small Higgs-boson masses, the cuts (fT2l) are almost as efficient 
as an experimentally not realizable cut on the pion energy in the x rest frame, as shown by 
the solid curves in Fig. [6ta). For heavy Higgs-bosons the discriminating power of the (p* 
distributions decreases, see Fig. Ob). This decrease can be avoided by an additional cut, as 



7 In fact, the searches 11231 12411 used pj > 20 GeV; but our tighter pj cut can, of course, always be applied in 



addition to the selected data sample. 



12 




Figure 7: Same as Fig. [6j including the contribution from the direct decay X — > % V. (a) m<j> = 
200 GeV; (b) m = 400 GeV. 

discussed in Sec. lIIICl 

In fact, the situation is improved by taking into account the contribution from the direct 
decay T~ — > %~ + v T which has maximal T-spin analyzing power. This decay channel is less 
strongly affected by the acceptance cuts (fT2l) because the n~ energy in the T~ rest frame 
is E K - — m z /2. This is shown by the ratio R of the contributions from the direct %~ and 
the p~ decays to the cross section pp 4> — >■ T~T + , given for m$> = 200 GeV and 
m<5 = 400 GeV in Table HB The numbers given in this table were computed at tree-level, 
but we expect them to not be strongly affected by radiative corrections to the 4> production 
amplitude. 

The (p* distributions, with the direct T~ — > %~ contribution included, are shown in Fig. |7J 
Comparing with Fig. [6] we see that the discriminating power has indeed improved, both for 
small and large Higgs-boson masses. For m# = 200 GeV, the cuts (fT2l) are in fact almost 
optimal, as one can see by comparing the dashed with the solid curves in Fig.|7J 
The reconstruction of the (p* distributions requires the determination of the (normalized) im- 
pact parameter vectors n_ and n + . One may ask whether a cut on their length would improve 
the sensitivity of the data selected in this way. We have therefore performed a simulation, 
along the lines outlined in where we require |n_| > 20 jlm for the displacement of the 
secondary vertex of T~ — > 7Z~, assuming an exponential decay of the x with a mean life-time 



ratio 


r,«l* = 200 
K-nocuts 


^/ttt = 200 

,r\ % - cuts 


7-,»i#=400 
I^nocuts 


n»J#=400 
P*~Mx- cuts 




0.31 


0.82 


0.31 


0.50 



Table II: Ratio of different final-state contributions to the cross section for pp — > <I> — > x x + with 
and without detector cuts as described in the text. 



13 



of 2.9 ■ 10 S. However, it turns out that the <p* distributions are only slightly affected - 
there is no gain in sensitivity. In addition, the cross section is reduced by almost a factor of 
2. Therefore, we refrain from this requirement in the following. 

C. Reconstruction of approximate T momenta 

From the discussion in the previous section we conclude that a more refined event selection 
is desirable for large Higgs-boson masses of the order of 400 GeV. In particular, additional 
cuts that remove events with low-energy pions (referred to the respective x rest frame), where 
b(E % ) is negative, would help to improve the discrimination of CP-even and CP-odd reso- 
nances. Knowledge of the T T 4-momenta would allow us to Lorentz-boost the measured 
pion momenta to the respective x rest frame, where the application of a cut on E K would 
be straightforward. An approximate reconstruction of the X momenta will be sufficient for 
this purpose, as long as it helps to enhance the difference of (p* distributions for scalar and 
pseudoscalar Higgs-bosons. In the following we describe an approach where we combine 
the information contained in the measured pion momenta in the pp laboratory frame and the 
experimentally known value of the Higgs-boson mass. 

We make the following approximations: i) In the laboratory frame, the x momenta k ± and % 
momenta p± are collinear, i.e., k ± = K±j) ± . ii) The measured missing transverse momentum 
■pmiss | s ass jg nec i to the sum of the differences of the transverse momenta of x ± and the 
charged prong a , i.e., P'™ ss = — p^ + — p£ . This is a very crude approximation for 
x — > p, a\, but it serves the goal formulated above. One can then write down eight equations 
for the eight unknown components of the x 4-momenta k ±fi : 



These equations can be solved analytically. One obtains an approximate Higgs momentum 
p^ 1 which can be used to boost the % momenta to the corresponding approximate Higgs- 
boson rest frame. We denote the resulting pion energies by . The distribution of 
is shown in Fig. H^a). The upper solid (black) curve shows the distribution for the p +a\ 
decays without cuts, while the dotted curves result from imposing the cuts (fT2l) . The lower 
solid and dotted (red) curves display the corresponding distribution for the direct x~ — > n~v 
decay. From this result and the analysis of the previous section we conclude that events with 
small pion energies in the x rest frames are correlated with events with small pion energies 
in the approximate Higgs-boson rest frame. This statement is supported by the distribution 
displayed in Fig.^b). This figure shows that, by imposing a cut on the distribution of 
E n - is shifted to larger values of the pion energy. 
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Figure 8: (a) Distribution of E~_ in the approximate Higgs rest frame. The solid curves correspond 
to the distribution without cuts; the dotted curves are obtained when applying the cuts Eq. (PT2l) . The 
vertical line indicates the additional cut E~_ = 60 GeV. (b) E n distribution for the combined decays 
T~ — > ,p~, with cuts (IT2T> . with and without an additional cut on E£_. 
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Figure 9: (a) The normalized <p* distributions (without the direct T~ — > %~ decay) with and without 
a cut on the reconstructed pion energy E^_ . The solid curves result from events with i? ~_ > 60 GeV. 
(b) Same as (a), but with the direct T~ — > n~ decay channel included. 



Encouraged by these observations we calculate, for heavy Higgs-bosons, the normalized <p* 
distributions by imposing the additional cut £21 > 60 GeV. The results are shown in Fig. [9j 
For comparison, the dashed curves result from applying only the cuts (fT2l) . The additional 
cut of > 60 GeV clearly leads to an increase in sensitivity. Including the x — > %v decay 
leads to a further improvement, see Fig. [9fb). In fact, the additional cut affects the direct 
decay channel only marginally. 

One should keep in mind that the additional cut on will reduce the size of the event 
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samples for the measurement of the <p* distributions. The sample based on T~ — > p , a\, and 
K~ decays will be reduced by about 18%. An analysis including a full detector simulation 
is required in order to optimize the selection cuts (PT2l) and the cut on ■ 
These conclusions will not be affected by analyzing different 4> production processes or by 
taking into account higher-order QCD corrections. For example, let us consider the produc- 
tion of a Higgs-boson with very high pj in the laboratory frame. Then the I and from 
T ± decays can pass the transverse momentum cuts ([111) and (fT2l) . even if the energies of the 
charged prongs in the respective T rest frames are small. As a result, contributions to the de- 
cay modes t — > /, a\ , p with b(E Jt j) > and b(E 7t j) < cancel and the discriminating power 
of the (p* distribution is reduced. However, when applying an additional cut in the approxi- 
mate Higgs rest frame as described above, the dangerous contributions will be rejected. This 
yields basically the same <p* distributions as before. In addition, for Higgs events with large 
Pt, much better methods of reconstructing the T rest frame can be applied, for instance the 
collinear approximation 1I27I1 or the method described in [28]. High pr particles/jets that are 
produced in association with a Higgs-boson can also be used to reconstruct an approximate 
Higgs rest frame. 



D. Combined leptonic and hadronic 1 -prong decays 

In this section we present results for the (p* distributions taking into account all 1 -prong 
decays ©, i.e., 

{ri'++x, l,l' = e,n, 
l~Z + + X and 7TI+ + X, (13) 
7Z-X++X. 

We refer to the different channels in (TT3T) by dilepton, lepton-pion, and two-pion final states. 
The cuts CCD and (PT21) are applied to the charged leptons and pions, / T and 7T T , respectively. 
For the dilepton final states the <p* distributions are presented in Fig. [TOfa) for two values 
of the Higgs-boson mass. The figure shows that the power of (p* to discriminate between 
a scalar and pseudoscalar Higgs-boson is, in these decay channels, almost independent of 
the mass of 4>. As discussed in Sec. IIII Al one would increase the sensitivity if one could 
reconstruct the T T rest frames and select an event sample with an additional cut on the lepton 
energies in these frames. However, Fig. [TOlb) shows that the enhancement would be rather 
modest. 

In Fig. fTTT a) and (b) the <p* distributions are presented for the lepton-pion and two-pion final 
states, respectively, for two different values of the Higgs-boson mass. For these final states, 
the discriminating power of the (p* distribution decreases with increasing Higgs-boson mass. 
We emphasize, however, that our evaluation is conservative in the sense that we applied only 
the acceptance cuts ([111) and (PT21) . As shown in Sec. IIII C[ a further cut on the charged- 
pion energy E£ in the approximate Higgs-boson rest frame would significantly enhance the 
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Figure 10: (a) The normalized <p* distributions for the dilepton final states, for m$ = 200 GeV 
(solid curves) and m$, = 400 GeV (dashed curves), (b) The solid curves are identical to (a) for 
m<& = 400 GeV; the dashed curves result from applying an additional cut E\ > ^ on the charged 
lepton energies in the x T rest frames. 
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Figure 11: (a) The normalized (p* distributions for the lepton-pion final states for m<$> = 200 GeV 
(solid curves) and m$ = 400 GeV (dashed curves), (b) The <p* distributions for the two-pion final 
states. 



discriminating power of (p* in the case of heavy Higgs-bosons or Higgs-bosons with large 
Pt- 

Notice that the (p* distributions for a scalar (pseudoscalar) Higgs-boson have opposite slopes 
for lepton-pion and two-pion final states. This is due to the fact that the signs of the leptonic 
and hadronic spin analyzer functions b(E{) and b(E n ) differ, both in the low-energy and 
high-energy part of the spectrum. Therefore, a very good experimental discrimination of 
leptons and pions will be crucial for this measurement. 
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mo [GeV] 


dilepton 


lepton-pion 


two-pion 


200 


380 


116 


18 


400 


600 


334 


207 



Table III: Event numbers needed to distinguish between a CP-even and CP-odd spin-zero state <I> with 
3 s.d. significance. 

It should be noticed that i) a Higgs-boson with scalar and pseudoscalar x- Yukawa couplings 
of equal strength (i.e., an ideal CP mixture) or ii) (nearly) mass-degenerate scalar and pseu- 
doscalar Higgs-bosons with equal production cross sections yield a _©* distribution which is 
flat, both for dilepton, lepton-pion, and two-pion final states (cf. yj]). In order to unravel 
these possibilities, one has to measure the distribution of the angle see below. The 
<p* distributions of mass-degenerate scalars and pseudoscalars with different reaction cross 
sections and of a CP mixture with \a T \ ^ |Z? T | have shapes which lie between the pure scalar 
and pseudoscalar cases and can also be disentangled with a joint measurement of the Ycp 
distribution. 

Next, we make a crude estimate of how many events are needed in order to distinguish 
between a scalar and a pseudoscalar Higgs-boson in the different decay channels (fT3l . We 
consider the asymmetry 

_ N(<p*>Jc/2)-N(q>*<Jc/2) 
r N((p* > 7t/2)+N((p* < jc/2) ' 

The asymmetries can be computed for the different final states from the distributions 
Figs. Qjla), CCDa), and (b), for a scalar and a pseudoscalar Higgs-boson. Assuming that 
systematic effects can be neglected, we estimate from these asymmetries the event numbers 
needed to distinguish a scalar from a pseudoscalar Higgs-boson with 3 standard deviation 
(s.d.) significance. These numbers are given in Table Hill 

One may ask how vulnerable the <p* distributions - and the distributions given in the 
next section - are with respect to uncertainties in the experimental determination of the 
<I> production/decay vertex and of the energies and momenta of the charged prongs. This 
question was investigated in |l] for the direct pion decay modes x~x + — > 7l^TZ + v x v x by a 
Monte Carlo simulation and it was found that these distributions retain their discriminating 
power when measurement errors are taken into account. One may assume that this result 
stays valid also for the larger class of 1 -prong decay modes considered in this paper. 
We close this section with a brief discussion of background reactions to the — > x~x + 
signal. These include QCD multijets, tt, single top, W + jets, Z/y* + jets, WW, WZ, and 
ZZ production. Among these, Z/y* — > x~x + (+ jets) constitutes an essentially irreducible 
background for a Higgs-boson with mass close to the Z mass. Most of this background can be 
distinguished from the signal by means of appropriate discriminating variables, in particular 
by reconstructing the T-pair invariant mass M xz using a likelihood technique, with which a 
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mass resolution of ~ 21% was achieved Q23D . For Higgs-bosons with masses m$ > 200 GeV, 
the Z* — > t~t + background can be suppressed by appropriate cuts on M xz , which works also 
for purely hadronic T x + decays I29I1 . 

If a Higgs-boson will be found with a mass not too far away from the Z mass one may, in the 
long run of the LHC, resort to the production of 4> by vector boson fusion, in order to study 
the above distributions. In vector boson fusion, 4> is produced in the central region, which 
provides a good veto against QCD background. More importantly, one has an additional 
signature from two well-separated forward jets, which gives a veto against Z* — > t~t + . Yet, 
in the SM and for large portions of the parameter spaces of models with an extended Higgs 
sector, gg — > <I> (and bb — > <& for large tan/3) is by far the dominant <I> production process. 
The modulus of the pseudorapidity, \r\ |, of a light Higgs-boson produced in these reactions 
is large, while its transverse momentum, generated by QCD radiation, is small on average. 
The XX pair from $ — > x~ x + is balanced in its total transverse momentum and its sum of 
the azimuthal angles. This provides a good discrimination against the QCD background, 
but not against Z — > t~t + . At this point one may exploit spin effects. The chiral-invariant 
Ztt and yxx couplings lead to characteristic XX spin correlations (see, for instance, Q30|-l33[]) 
which differ from those that result from the decay of a spin-zero resonance whose fermion 
couplings are chirality-flipping. For instance, for the decays T~T + — > 7T~7T + vv, this has 
the following consequence. If the n~iz + result from Z boson decay, the number of 
events with E K - and E K + both large or both small (in the t~t + ZMF) is much larger than 
the number of events with E K - large (small) and E K + small (large); while for <I> — > t~t + 
just the opposite is the case. These %% energy distributions may be used to discriminate 
the signal from the irreducible background in the case of a light <I>. In addition, also the 
7t~K + invariant mass distribution shows some difference between events from 4> — > t~t + 
and Z — > t~t + I32Q. Rather than trying to discriminate against the irreducible background, 
an alternative strategy might be, for a light <I>, to take into account the Z — > x~ x + events both 
in the measurement and in the Monte Carlo modeling of the distributions © and (flOl) . This 
requires a detailed study which is, however, beyond the scope of this paper. 



E. Higgs-sector CP violation 

Besides <p*, a further important observable in this context is the angle Xf/^p defined in (flOl) . 
It is the appropriate variable to check whether or not a spin-zero resonance <I> has couplings 
to both scalar and pseudoscalar x lepton currents. A nontrivial \\Tq P distribution, respectively 
a nonzero asymmetry associated with this distribution would be evidence for CP violation 
in the "Higgs sector" (which is different from Kobayashi-Maskawa CP violation). Such a 
discovery would have enormous consequences, in particular for baryogenesis scenarios (see, 



e.g., the reviews 041 I35IR 



We assume here that 4> is an ideal mixture of a CP-even and CP-odd spin-zero state, with 
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Figure 12: The i//£ p distributions for (a) dilepton, (b) lepton-pion, and (c) two-pion final states. The 
different scenarios are explained in the text. 

reduced Yukawa couplings a T = —b x to X leptons 8 . For definiteness, we take a x = —b x = 1. 
We call this the CPmix scenario for short and consider it for a Higgs-boson with mass m<j> = 
200 and 400 GeV. For comparison we consider also three scenarios where CP is conserved: 
i) a pure scalar H, ii) a pure pseudoscalar A, and iii) the case of a (nearly) mass -degenerate 
scalar H and pseudoscalar A with approximately the same production cross section and decay 
rate into x leptons. 

The distribution of Xjr^p is displayed in Figs. [T27 a) - (c) for dilepton, lepton-pion, and two- 
pion final states, respectively. The figures show that the variable ^ p efficiently distinguishes 
between CP conservation and violation - for the CP-conserving scenarios i) - iii) above, the 
distribution is fiat. For a CP-mixed state, the slope of the i/z^p distribution for lepton-pion 

8 Suffice it to mention that for the normalized y/^p distribution only the relative magnitude and phase of a T 
and b T matter, while the magnitudes of these couplings determine the decay rate of <J> — > XX. 
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[GeV] 


dilepton 


lepton-pion 


two-pion 


200 


1540 


514 


116 


400 


2400 


1390 


726 



Table IV: Event numbers needed to find evidence with 3 s.d. significance that <I> is an ideal CP mixture. 



final states is opposite in sign to the slope for dilepton and for two-pion final states, for 
reasons mentioned above. In analogy to (fT4l) one may consider the asymmetry 

_ #(yfo> K/2) ft/2) 
N(\j/* p > ft/2) +N{\j/* p < ft/2)' 1 ' 

With the values of Ay* obtained from the distributions Figs. \T2\ a.) - (c) one gets the esti- 
mates of the event numbers, given in Table HVl that are needed to find evidence with 3 s.d. 
significance that $ is an ideal CP mixture. 

One may enhance the discriminating power of the t//^ distribution by constructing an ap- 
proximate Higgs-boson rest frame as outlined in Sec. MI Cl and impose additional cuts on the 
energies of the charged-pions in this frame. In analogy to the above results for the (p* distri- 
bution, we find that the improvement as compared to the results shown in Figs.[T27b'). (c) is 
small for Higgs-boson masses below 200 GeV, whereas it becomes significant for heavy CP 
mixtures with m& ~ 400 GeV. 



IV. CONCLUSIONS 

We have shown that the CP quantum numbers of a Higgs-boson 4>, produced at the LHC, 
can be determined with the observables © and (flOl) in the T-decay mode ^> — > T + T~, using 
all major subsequent 1 -prong x decays. The selection cuts that we applied in our analysis 
to the dilepton, lepton-pion, and two-pion finals states significantly enhance the discriminat- 
ing power of these observables for the "non direct charged-pion decay modes". Therefore, 
practically all x decay modes can be used for pinning down the CP properties of because 



the three-prong x decays can also be employed for this purpose II 1811 . Depending on the 4>- 
production cross sections, i.e., on its mass and couplings, it should be feasible to collect the 
event numbers (estimated in Tables Hill and HVl) that are required for statistically significant 
CP measurements after several years of high-luminosity runs at the LHC. 



APPENDIX 



In this appendix we collect, for the convenience of the reader, some results on x decays which 



are relevant for the calculations described above. (For a review, see |36| 
ratios of the 1 -prong T-decay modes, given in Table El are taken from Q37 



.) The branching 
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decay mode 


x ± ^n ± 


X ± — > p ± — > Tt^K 


-»• af -»• w ± 2«° 


x ± -^e ± ,n ± 


BRpdg [%] 


10.91 


25.51 


9.3 


35.2 



Table V: Branching ratios for the major 1 -prong T-decay modes 13711 . 



Next we list the spectral functions n(E a ) and b(E a ) of the energy-angular distributions © of 
polarized x^ decays to a T . The functions given below apply to both T~ and x + - but notice 
the sign change in front of b(E a ) in ©. Furthermore, our convention for the distribution © 
is such that we differentiate with respect to the energy E a of the charged prong. Therefore, 
the functions n[E a ) are dimensionful while the functions b(E a ) are dimensionless. 



The decay x T ->■ 7i T + v T 



In the 2-body decay x — y % -\- Vx the energy E 1 ^ in the x rest frame is fixed and the functions 
n n (E n ) and b n (E n ) are given by Q38[]: 



2m T 



M £ w) = 1. 



(16) 



The decay x 1 



The differential rate of the decay of polarized x leptons to a charged pion via a p -meson was 
calculated in 113 811 . With x = 4E n /m t , where E n denotes the energy of the charged-pion in 
the x rest frame, the spectral functions are given by 

6 (x-r-\) 2 + (l-r)(r-p) 



n p (E K ) 
bp (Ex) 



m T (l-r) 2 (l + 2r)(l-p/r) 3 /2' 

x(x — r— l) 2 +x(3 — r){r — p) —4(r — p) 



x 2 -4p((x-r- 1) 2 + (1 -r)(r-p)) 



with p = 4m 2 /m 2 and r — m 2 /m 2 . These functions are plotted in Fig. Ufa). The kinematic 
range of E n is 



* i + ,-(i-,)./i-£ )<e k <^ l l+r+(l - r) M-P 



(17) 



The decay x T ->■ a f ->■ 7r ZF 27r° v T 



The differential rate of the 1 -prong decay of polarized T leptons to a charged-pion via a 
ai-meson was calculated in 113911 . The corresponding functions n ai (E K ) and bn x (Ex) are 
complicated and were fitted to the numerical results shown in Fig. 6-4 of Ref. [13911 . With 

2m x (E % -m n ) 



x = 



m 2 — 3m 2 — 2m x m % 



(18) 
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where E K is the energy of the charged-pion in the x rest frame, we obtain 

n a AE„) = - - - x — (0.0112624-2.15495*+ 165.368* 2 

m\ — 3m^. — 2m x m K \ 

-997.586/ + 2818.75/ -4527.77/ (19) 
+4250.43/ - 2182.33/ + 475.283/ 

b ai (E n ) = -5. 28726Vx + 9.38612jc- 1.26356/ (20) 
-18.9094/ + 36.0517/ - 19.41 13/ . 

The plots of n a] (E n ) and b ai (E n ) are shown in Fig. HJb). The kinematical range of the 
charged-pion energy in the x rest frame is 

ml — 3ml 

m n <E n <^- n -. (21) 

2m T 

The decay x T ->■ l T ViV t 

For the leptonic decays x ± — > / ± V/V T the mass of the final state lepton, e or u, can be ne- 



glected. Using x = 2Ei/m T , where Ei is defined in the x rest frame, one has 113 811 



= — / (3 - 2x) , btiEt) = (22) 



with < Ei < m t /2. 
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